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Abstract 

We illustrate in a simple setting the instantaneous shock tracking approach to sta- 
bility of viscous conservation laws introduced by Howard, Mascia, and Zumbrun. This 
involves a choice of the definition of instanteous location of a viscous shock- we show 
that this choice is time-asymptotically equivalent both to the natural choice of least- 
squares fit pointed out by Goodman and to a simple phase condition used by Gues, 
Metivier, Williams, and Zumbrun in other contexts. More generally, we show that it is 
asymptotically equivalent to any location defined by a localized projection. 



1 Introduction 

In this note, we illustrate in the simple and concrete setting of Burgers equation the ar- 
gument for nonlinear stability of viscous shock waves developed for general systems of 
conservation laws in [ZH IMaZ21 IMaZ31 IMaZ4] , based on instantaneous tracking of the lo- 
cation of the perturbed viscous shock wave. The advantage of Burgers equation is that 
the linearized equations may be solved explicitly by a linearized Hopf-Cole transformation, 
thus isolating the nonlinear issues we wish to discuss. This same example was given in |Z1] : 
here we expand a bit the surrounding discussion, reexamining the question of what is a 
reasonable or natural definition of the instantaneous location of a perturbed viscous shock 
wave and adding a discussion of the small- amplitude limit. 

Using the purely operational but analytically tractable definition of [Zl^ IMaZ2l IMaZ4] 
as a tool for comparison, we show that any definition based on localized projection is 
time-asymptotically equivalent to any other and to the definition of [Zl^ IMaZ2t IMaZ4j : 
see Appendix [HI and especially Remarks IB.3f[B74l Moreover, any of these may be used as 
the basis of a nonlinear stability argument. This includes in particular both the natural 
definition by least squares fit pointed out early on by Goodman [G] and, in the limit of 
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infinite localization (to a single point), a very simple definition based on a phase condition, 
introduced by Gues, Metivier, Williams, and Zumbrun [GMWZ] . 

Our analysis is intended for the nonspecialist. It is brief and self-contained except for 
standard linear and short-time parabolic existence theory. Though we restrict here for sim- 
plicity to the scalar Burgers case, our arguments and conclusions extend in straightforward 
fashion to the general system case \Z1\ IMaZ21 IMaZ4| . once there are established the req- 
uisite bounds on the linearized solution operator. This separate, and in general difficult, 
problem has been treated in |ZH1 IMaZ31 IZ2j : see Remark 13.41 and the discussion of Section 
m Our purpose here is, rather, to isolate the issues connected with viscous shock-tracking 
and the nonlinear iteration argument by restricting to a case where the linearized bounds 
are available by exact solution formula. 

1.1 Problem and equations 

Consider the scalar viscous conservation law 

(1-1) Ut + f{u)x = Uxx, 

u = u{x, t) eR, X £R,t£ M+, with 

(1.2) f{u) = nV2. 

Eq. (II. ip serves as a simple model for gas dynamics, traffic flow, or shallow-water waves, 
where u represents the density of some conserved quantity and / its flux through a fixed 
point X. With the choice of flux (jl.2p . (jl.ip becomes Burgers equation, the prototypical 
example of a scalar viscous conservation law. Behavior for other (convex) fluxes is qualita- 
tively similar. 

We investigate the question of nonlinear stability of solutions u = u, that is, whether a 
perturbation u = u + u remains close to u in some norm for initial perturbations u\t=o = 
{u — u)\t=Q sufficiently small in some (possibly different) norm: more specifically, nonlinear 
asymptotic stability, that is, whether u both remains near and converges to u as t ^ +oo 
for initial perturbations sufficiently small. Since the equation (II. ip is translation-invariant, 
we must when relevant (specifically, when translates of u are not equal to u) adjust the 
second notion to that of nonlinear asymptotic orbital stability, defined as nonlinear stability 
together with convergence as t — > -|-oo to the set of translates of u, as discussed further 
below. 

1.2 Constant and traveling- wave solutions 

An obvious class of solutions of (jl.ip are the set of constant solutions u = a, a G M. A 
second class of solutions are viscous shock waves, or smooth traveling-wave solutions 

(1.3) u{x , t) = u{x — st) , lim u{x) = u± 
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s constant, connecting constant endstates u±. If s = 0, they are equilibria, or stationary 
waves of the associated evolution equation (jl.ip . A traveling wave may always be converted 
to a standing wave by the change of coordinates x ^ x — st to a frame moving with the 
same speed s. 

Observing that dtu{x — st) = —su', dxu{x — st) = -u', and d'^u{x — st) = u" , we obtain for 
a solution (jl.3p the profile equation —su' + /{u)' = u" . Integrating from — cxd to x reduces 
this to a first-order equation 

(1.4) u' = {f{u)-su)-{f{u^)-su.). 

For definiteness taking s = 0, u_ = 1, we obtain u' = (1/2)(1 — v?), which has the explicit 
solution 

(1.5) m(x) = -tanh(x/2) 

connecting endstates u± = ^1. This is the unique solution up to translation in x connecting 
that particular pair of endstates. Other endstates and speeds also lead to tanh profiles, as 
may be seen by invariances of Burgers equation; thus, we may without loss of generality 
restrict to this specific case. 

2 Stability of constant solutions 

To indicate the basic approach, let us first consider stability of a constant solution 

(2.1) u = a, aeR 

of (II. ip . Letting n be a second solution of (jl.ip . and defining perturbation u := u — u, we 
obtain after a brief computation the perturbation equation 

(2.2) ut-Lu = N{u)^, 

where Lu := the linearization of Uxx ~ f{u)x about solution u — a. and 

N{u) := —v} 12 is a quadratic order remainder. 

2.1 Linear solution operator 

The homogeneous linearized equations vt — Lv = may be recognized as a convected heat 
equation 

(2.3) vt + avcc = vxx, v\t=o = f. 
This admits an exact solution 

/ + 00 
Gix,t;y)f{y)dy, 
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where 



(2.5) 



_ \x-y-at\ 

G{x,t;y) := e^%(x) = ■ 



is the Green function for p.3p . a convected heat-kernel. This yields in particular a unique 
classical solution v G C°(t > 0; LP{x)) n C^it > 0, x) for each / e L^. 
Easy scaling arguments yields, for 1 < p < oo, 



(2.6) 



|G(-,t;y)|LP(,) = |G(x,t;-)|Lp(,) =Cpt-^(i-i/^'), 



for some constants Cp, C'p > 0. From (|2.6p . we readily obtain the following linearized 
estimates (standard heat kernel bounds). 



Lemma 2.1. For some C > 0, all t > 0, 

" + CXD 



(2.7) 



oo 

+ 00 



G{x,t;y)f{y)dy < Ct-^^^-^/^') 
Gy{x,t;y)f{y)dy\^^ < Ct-^d-Vp)"! |/|^,^^^. 



Proof. Applying the Triangle inequality together with (j2.6p . we obtain 



+ 00 



G{x,t;y)f{y)dy 



LP{x) 



< 



\Gi;t;y\Lr>i.)\f{y)\dy = 



□ 



The proof of the second inequality is similar. 

Lemma 2.2. For some C > 0, all t > 0, 
(2.8) 

/+00 r+oo 
G{x,t;y)f{y)dy ^ < / Gy{x,t-y)f{y)dy ^ < 

Proof. Noting that G{x, t; y) = G{x — y, t; 0), so that (j2.4p is a convolution, we may rewrite 
I-^ G{x,t;y)f{y)dy as J^^ G{z,t;0)f{x - z)dz with z := x -y. Applying the Triangle 
inequality and (|2.6p . we obtain 



G{x,t-y)f{y)dy < \G{z,t;0\\f\LPdz = Gi\f\LP. 

LP(x) 



The proof of the second inequality is similar. 



□ 
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2.2 Integral representation 

From the homogeneous hnearized solution formula (|2.4|) . we obtain by variation of con- 
stants/Duhamel's formula a solution for the inhomogeneous linearized equations 

(2.9) vt- Lv = g, v\t=o = f 

oiv = e^7 + /o e^^^-'^g{s)ds, or 

/+00 />+oo 
G{x,t;y)fiy)dy+ / / Gix,t - s;y)giy, s)dy ds, 
"Oo Jo J —oo 



.P 



yielding a unique C°(t > 0;LP(x)) n C^(t > Q;x) solution u for / G and 5 G 
2.3 Nonlinear iteration 

Returning now to the nonlinear problem (|2.2p . we have, setting = N{u)x in (|2.9p . the 
representation u{x,t) = f^^ G{x,t;y)uo{y)dy + /q - s;y)N{u{y, s))ydy ds, or, 
integrating the last term by parts, 

G{x,t;y)uoiy)dy - / Gy{x,t - s;y)N{u{y, s))dy ds, 
-00 Jo J —00 

valid so long as the solution u exists and remains sufficiently smooth that (j2.1ip gives the 
unique solution to the associated inhomogeneous problem, in particular for uq € H L°° 
and u in C°{t > 0; n L°°) H C^it >0,x), any p > 1. 

On the other hand, standard short-time existence theory (proved, e.g., by contraction- 
mapping using a similar representation with shifted initial time) yields existence of aC^{t > 
0; LP n L°°(a;)) nC^{t> 0, x) solution so long as |u|iPnL°° remains bounded. 

Define now 

(2.12) C{t):= sup |n|Lp(s)(l + f)5(i-i/p). 

0<s<t, l<p<oo 

Lemma 2.3. For all t >0 for which (^(t) is finite, some C > 0, and Eq := |?io|Linl,°°7 

(2.13) at)<G{Eo + atf). 

Proof. Noting, by quadratic dependence N{u) = 0{\u\'^) and the definition ()2.12p of that 

^2^^) \N{u)\Li<G\u\l,<C{tf{l + t)-h 

\N{u)\lp < C\u\lp\u\l^ < C{tf{l + t)-i(i-i/p)-|, 



5 



we obtain, applying Lemmas l2.7H2.8l to representation (j2.1ip . the estimate 

" + CXD 

G{x,t;y)uo{y)dy 



< 



t/2 p+oo 



(2.15) 



Gy{x, t - s; y)N{u{y, s))dy ds 

/•t f+oo 

/ / Gy{x,t - s;y)N{u{y,s))dyds 

J t/2 J-oo 

ft/2 

< C(l + t)-i(i-i/p)Eo + CCitf / {t- s)-|(i-i/p)-i/2(i + s)-Us 

Jo 



LP{x) 
LP{x) 



+ CC{tf {t- + s)~^(l-l/P)-3cis 

Jt/2 



Rearranging, we obtain (12.31 



□ 



Corollary 2.4 (Stability of constant solutions). Constant solutions u = a are nonlinearly 
stable in n L°° and nonlinearly asympotically stable in LP , p > \, with respect to initial 
perturbations uq that are sufficiently small in n -L°°. More precisely, for some C > 0, 

(2.16) \u-u\Lv{t) < 0(1 + t)-'2(^-^/p)\u - u\LinL^\t=o 

for allt>0, 1 < p < oo, for solutions u of (II. ip with \u — u\L^r\L°°\t=o sufficiently small. 
Proof ("Continuous induction") By Lemma ({t) < C{Eo + C{ty) for 

(2.17) £^0 := \u - n|iinLoo|t=o- 

Taking Eq < -^tj, we have therefore that C,{t) < 2CEq whenever C,{t) < 2CEq, and so the 
set of t > for which C{t) < 2CEo is equal to the set of t > for which C{t) < 2CEo. 
Recalling, by the cited standard short-time existence theory, that C is continous wherever 
it is finite, we find, therefore, that the set of t > for which C,[t) < 2CEq is both open and 
closed. Taking without loss of generality C > 1/2, so that t = is contained in this set, we 
have that the set is nonempty. It follows that C(t) < 2CEq for all t > 0, yielding (12.16P by 
definitions (l2T2]l and (|2T7|) . □ 



Remark 2.5. The rate of decay (j2.16p is that of a heat kernel- that is, the mechanism for 
stability is diffusive only. 

3 Stability of viscous shock solutions 

We turn now to the stability of viscous shock solutions of (jl.ip . without loss of generality, 
restricting to the case 

u{x) = — tanh(x/2) 
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described in (jl.Sp . Letting u as before be a second solution of (jl.ip . define the perturbation 

(3.1) u{x , t) := u{x + a{t) , t) — u{x) 

as the difference between a translate ofu and the background wave it, where the translation 
a{t) is to be determined later. 

This yields after a brief computation the perturbation equation 

(3.2) ut - Lu = N{u)::, + a{t){ux + Ux), 

where Lu := Uxx — {0'{x)u)x is the linearization oiuxx—f{u)x about solution u = — tanh(x/2), 
a{x) := df{u){x) = u{x), and N(u) := —u'^/2 is the same quadratic order remainder as in 
the constant-coefficient case. 



3.1 Linear solution operator/decomposition of the Green function 

The homogeneous linearized equation 

(3.3) vt- Lv = vt + ia{x)v)x - Vxx = 0, v\t=o = f 

can be solved explicitly by linearized Hopf-Cole transformation [Si IN ] [Z3 l IGSZ| . to give an 
exact solution formula 

/ + 00 
Gix,t;y)fiy)dy, 
-oo 



where 



(3.5) 



G{x, t; y) := e^''5y{x) = u'{x) (errfn( ^ J_ - errfn( ^ J^^ ' 



e^+e"^^ \/47rt \e^+e~^' \/4vrt 

is the Green function for ()3.3p and errfn(2:) := ^ e'^^d^. 
Following the approach of \Z1 \ IMaZ"2l IMaZ4] . decompose now 

(3.6) G(x, t; y) := E{x, t; y) + S{x, t; y) + R{x, t; y), 
where 

(3.7) E{x, t; y) := n'(x)e(y, t), e(2/, t) := (^) (errfn(^^) - errfn(^^)) , 



(3.8) S{x, t; y) := (( / ' . ) " + ( . . ) " " ) , 

^^62+6 2/ v4'/rt ^62+6 2/ v4vrt ^ 
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and 



(3.9) ^ ^ 

— ■u'(x) ( - ) ( errfn(^^^ — ^= — ) — errfn( ^ 



i2(x, i; y) := ti'(x) (-) (errfn( ^= — ) — errfn(- ^ 



2/V ' V4t V4t 

Here, the "excited term" E represents the nondecaying part of the hnearized solution v, 
involving the zero-eigefunction Lu' = associated with instantaneous translation of the 
background wave, the "scattering term" S comprises Gaussian signals convected along hy- 
perbolic characteristics, and the "remainder term" R a faster-decaying residual. 
A straightforward calculation gives 



(3.10) 



\R{x,t;y\<C\x\\u'ix)\ I ( — ^|^| J^^"^"^^ 



^ 

<Ce-^lV r A + "~ A )d9, 
Jo ^ VAirt V4vrt 



9 > 0, showing that R, as the product of an exponentially decaying term and the sum of 
convected Gaussians, is indeed faster-decaying than either E or S. 

Lemma 3.1. For some C > 0, 9 > 0, allt > 0, 



(3.11) 



\R{x,t;y\ < Ce-'\^\l^{ ^ + ' ) + Ce~'^\^~y\+'\ 



\Ryix,t;y\ < Ce-^N/C(^- + ^- + 



t / Vt 



Proof. Applying the Cauchy-Schwarz inequality in the argument of the exponential, we 
find readily that 

(3.12) e ^l^l/^fe 4t +e ~t ]<e ^1^1/'-' (e ici +e 4ct 

for |x| < Mt and C > sufficiently large, hence (IS.lOp implies (I3.1ip (i). For \y\ » |x| -|- |t, 
(j3.12p holds trivially, likewise giving (j3.1ip (i). In both of these cases, the lefthand side is 
bounded by the first, Gaussian, term alone on the righthand side. In the remaining case 
|x| >> t and \y\ < M\x\, we have for C > sufficiently large that e"^'^' < e~*^''/'-^)'^l^~2/l+t)^ 
from which we find directly from ()3.9p that the lefthand side of (|3.1ip (i) is bounded by the 
final term on the righthand side. 

Similar computations yield (j3.1ip fiil. □ 

Remark 3.2. The excited term E converges as t ^ -|-oo to u'{x) times 

;>+oo r+oo 

-ai+oo):= e(y,+oo)/(y)dy = (1/2) / f{y)dy, 
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the time-asymptotic state of the hnearized equations (j3.3p determined by conservation of 
mass (equals total integral v{x,t)dx). Note that u'{x) corresponds to infinitesimal 
translation of the background wave u{x), hence a linear time-asymptotic state —(tu'{x) 
corresponds roughly to a steady-state perturbation u{x — a) — u{x) consisting of a shift, 
or translation, a of the background wave. The term a{t) := — e{y,t)f{y)dy thus 
measures, at a linearized level, the shift in location of the shock at time t, or "instantaneous 
shock shift" . This refines the picture of behavior given by the time-asymptotic shock shift 
cj(+oo). 



Proposition 3.3. The Green function G decomposes as G 
where, for some C > 0, all t > 0, 



E + G, E = u'{x)e{y,t), 



(3.13) 
(3.14) 



G{x,t;y)f{y)dy 



— oo 

+ 00 



Gy{x,t;y)f{y)dy 



LP(x) 



LP{x) 



<Gt''2(^''M\f\^^, 
<Ct-i(i-iM-||/|^,. 



f + 00 



G{x,t;y)f{y)dy 



LP{x) 



<G\f\LP, 



+ 00 



Gy{x,t;y)f{y)dy 



LP(x) 



<Gt-2\f\LV^ 



(3.15) 



etiy,t)fiy)dy 



<Gt-2\f\Li 



ety{y,t)f{y)dy <Gt-'\f\Li 



(3.16) 

and 

(3.17) 



+ 00 



etiy,t)f{y)dy <G\f\ 



e{y,t)f{y)dy < 



eytiy,t)fiy)dy < Ct '^\f\L--] 



{y,t)f{y)dy <Ct-i/VlLi- 



Proof. Defining G := R + S,we have the decomposition G = E + G. By (j3.8p and estimate 
()3.1ip . G and Gy obey essentially the same bounds as G and Gy in the constant-coefficient 
case ()2.5p . up to a harmless exponential error (the final terms on the righthand sides of 
(j3.1ip l. Thus, bounds (j3.13p and (j3.14p follow by the same argument used to prove (j2.7p 
and ()2.8p . By (j3.7p . |et| and \eyt\ satisfy essentially the same bounds as sup^ IGj and 
sup^. \ Gy\, hence (|3.15p and (|3.16p follow again from this same argument in case p = oo, 
which amounts to Holder's inequality together with bounds on e and derivatives (see 
Lemma IC.H Appendix IC.ll for a careful derivation of these bounds). Finally, (j3.17p 
follows by |e| < C, \ey\ < Gt~^^'^ using the triangle inequality. □ 
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Remark 3.4. The apparently special Proposition 13.31 in fact holds for viscous shock waves 
of general strictly parabolic systems provided that the shock satisfies a generalized spec- 
tral stability, i.e., Evans function, condition \Z1\ IZ41 IMaZ3j . Indeed, there is a parallel 
decomposition of the Green function as the sum of terms E, S, and R with pointwise de- 
scriptions generalizing those of (|3.7|) . (|3.8|) . (|3.11|) . Similar bounds hold for Evans stable 
shocks of general hyperbolic-parabolic systems [MaZSt [Z4] . Scalar shock waves are always 
spectrally stable, by the maximum principle; hence, the stability condition does not make 
itself apparent for Burgers equation. 

In the derivation of bounds by inverse Laplace transform estimates, the terms E and 
S arise in a very natural way as leading terms of a low-frequency "scattering" expansion 
|MaZ3[ IZ2] of the resolvent kernel about frequency A = 0, without the need to re- arrange 
terms as done here in the Burgers case. See Section 2, [BeSZj . for a particulary clear 
discussion of the method from more general point of view. Indeed, the decomposition of 
G into E and G was suggested from the inverse Laplace transform point of view [ZHl IZl^ 
IMaZ3| . Here, for pedagogical purposes, we have imposed this structure by force on the 
explicit Green function given by Hopf-Cole transformation in order to demonstrate clearly 
the approach. 

3.2 Integral representation/a-evolution scheme 

Recalling that u'{x) is a stationary solution of the linearized equations ut = Lu, so that 
Lux = 0, or 

roc 

G{x,t;y)ux{y)dy = e^^Ux{x) = Ux{x), 



we have, applying Duhamel's principle to (j3.2 



u{x,t)=l G{x,t;y)uo{y)dy 

(3.18) 



— oo 
t />oo 



Gy{x, t — s; y){N{u) + au){y, s) dy ds + a{t)u{x) 
Defining a implicitly as 

/oo 
e{y,t)uQ{y)dy 

+ / / ey{y,t - s){N{u) + au){y,s)dyds, 



J-oo 



following [ZHl EH IMaZ2l IMaZ3j . where e is defined as in (f37fl) . and substituting in (f3?T8]) 
the decomposition G = u'{x)e -|- G of Proposition 13.31 we obtain the integral representation 



u{x,t) = I G{x,t;y)uo{y)dy 

(3.20) 



oo 

t fOO 



y 

J-oo 



Gy{x, t — s; y){N{u) + au){y, s)dy ds, 
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and, differentiating (j3.19p witli respect to t, and observing tliat ey{y,s) — > as s — > 0, as 
the difference of approaching heat kernels, 



a{t) = - J et{y,t)uo{y)dy 

(3.21) r+oo 

+ / eyt{y,t - s){N{u) + au){y,s)dyds. 



J-oo 

Equations (j3.20p . (j3.2ip together form a complete system in the variables («,«), from 
the solution of which we may afterward recover the shift a via (|3.19p . From the original 
differential equation (j3.2p together with (|3.2ip . we readily obtain short-time existence and 
continuity with respect to t of solutions {u, a) £ n L°° x M by a standard contraction- 
mapping argument 1^ 

Remark 3.5. Here, the key step in deriving (j3.20p is to observe that the contribution in the 
righthand side of (j3.18p coming from terms involving u' {x)e{y,t) is, under definition (j3.19p . 
exactly —u'{x)a{t), so cancels the final term. That is, we have defined the instantaneous 
translation a{t) from considerations of technical convenience so as to cancel all nondecaying 
terms in (|3.18p . Note that a{t) agrees to linear order with the prescription a{t) in Remark 
of the instantaneous shock shift for the linearized equations. 



3.3 Nonlinear iteration 

Associated with the solution {u,a) of integral system ()3.20p - (l3.2ip . define 

(3.22) Q{t):= sup {\u\lp{s){1 + tp^~^/'''^ + \a{s)\{l + sf'^) . 

0<s<t, l<p<oo 

Lemma 3.6. For all t > for which C(t) is finite, some C > 0, and Eq := |no|iinL°°; 

(3.23) at)<C{EQ + atf). 

Proof. With the established bounds on G and e, the proof of (|3.23p is almost identical to that 
of (12.13P in the constant-coefficient case. Noting, by quadratic dependence N{u) = 0(|np) 
and the definition (j2.12p of Q, that 



\N{u) + au\Li <C\u\Li{\u\L^ + \a\) <Q{tY{l + t 



(3.24) 

\N{u)+au\Lv < C\u\lv{\u\l^ +\a\) < C(t)^(l + t)" 2 (i-Vp)- 5 ^ 



^ Specifically, for initial time T > 0, and t > T, split the expression (I3.21|l for a{t) into the sum of a 
bounded "known" term — f^^et{y,t)uo{y) dy+ eyt{y,t~ s){N{u)+au)(y, s) dy ds and an "unknown 

term" ^yt{y, i — s){N{u) + au){y, s) dy ds that is contractive for (u, a) bounded and \t ~ T\ << 1. 

The u-equation (|3.2|l may be treated in standard fashion, treating the righthand side as a forcing term and 
expressing u as an integral on [T, t], again contractive for |i — T] << 1. 
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we obtain, similarly as in (j2.15p . applying Lemmas l3.13fl3TT^ to representation (j3.20p . 

< C{1 + + cat? f'\t - ,)-|(l-l/p)-l/2(i + 

JO 

(^•2^) +catf f + 

Jt/2 

< C{EQ + Q{tf){l+t)-^'^^-^'P\ 
Similarly, by (l3l^ and (|3l6]) . 

\6i{t)\ < C(l + t)-2£;o + CCitf / (t - s)-\l + s)-2ds 



(3.26) , ^^,,,2 / , 



+ CC(t)' / {t- s)-Hl + s)-'ds 

Jt/2 

<c{Eo+atf){i+t)-h. 

Combining and rearranging (|3.25p - (|3.28p . we obtain (|2.3p . □ 

Corollary 3.7 (Stability of shock solutions). Viscous shock solutions u{x) of (jl.ip are 

nonlinearly stable in H and nonlinearly orbitally asympotically stable in LP , p > 1, 
with respect to initial perturbations uq that are sufficiently small in L^r\L°° . More precisely, 
for some C > and a G W^'°°{t), 



\u - u{- - a)\Lp{t) < C(l + t)-5(i-i/p)| 



^ ~ '"lLinL°°|i=0) 



(3 27) ^ C'(l + i)"2|u-u|iinLoo|t=o, 

\a{t)\ < C\u - u\L^nL^\t=o, 
\^ - n|iinL«=(t) < C|n - u\L^nL^\t=o, 

for allt>0, 1 < p < oo, for solutions u of (ll.ip with \u — u\L^nL°°\t=o sufficiently small. 

Proof. The first two inequalities follow by a proof identical to that of Proposition 12.41 in the 
constant-coefficient case, using (13. 6p and continuity of ^ wherever is finite, a consequence 
of short-time existence theory, to obtain ({t) < 2CEq, for Eq := \u — n|^i|-|^cx; |t=o < "Ho 
sufficiently small. This yields the first two bounds by definition of C. The third then follows 
using (f3T7p . by 



t/2 



\a{t)\ < CEo + CCity / {t - s)-2{l + s)-2ds 
Jo 

^^•^^^ + CCitf f {t - S)-5(1 + s)-5ds 

Jt/2 

<c{Eo+at?)- 
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Finally, we note that 

u{x, t) — u{x) = u{x — a{t),t) + {u{x) — u{x — a(t)), 

so that \u{-,t) — u\ is controlled by the sum of \u\ and \u{x) — u{x — ct{t))\ ~ a{t)\u' {x)\, 
hence, by our estimates, remains < CEq for all t > 0, for £^0 sufficiently small. This verifies 
the fourth inequality, yielding nonlinear stability and completing the result. □ 

Remark 3.8. In the semilinear case considered here. Corollary 13.71 could be proved in 
more straightforward fashion by a contraction mapping argument applied directly to the 
system (I3.20p -( f"3.2ip . bypassing the continuous induction argument above. However, in 
more delicate situations such as the quasilinear parabolic or hyperbolic-parabolic case, it is 
advantageous for reasons of regularity to separate the issues of short-time existence/well- 
posedness and long-time bounds, as we have done here; see |MaZ21 IMaZ41 IZ41 IRZj for 
further discussion. 

Remark 3.9. Again, the rate of decay (lA.lSp is that of a heat kernel- that is, the mecha- 
nism for stability is diffusive only, and not involving compressivity of the shock. This rate is 
in fact sharp, as may be seen intuitively by considering a compactly supported perturbation 
supported arbitrarily far from the shock location x = 0. Far from the shock, the background 
solution u is approximately constant, and so behavior is like that of a perturbation of a 
constant solution as studied in Section [2l But, this is readily seen to decay like a heat 
kernel, giving the stated rate (jA.lSp . 

3.4 Postscript: phase-asymptotic vs. asymptotic orbital stability 

A stronger condition that nonlinear orbital stability, proved above, is nonlinear phase- 
asymptotic orbital stability, in which a perturbed solution u is required to approach not 
only the set of translates of u, but a specific translate of u. In the language of Corollary 
13.71 this amounts to the requirement that a have a limit a{t) a{+oo) as t ^ +00. 

We do not establish this property in Corollary 13.71 nor is it estabhshed in [ZH IMaZ 2t 
IMaZ4| . Indeed, for the general class of perturbations considered here (and in [Zll, iMaZ^ 
IMaZ4j ). a(t) if it converges to a limit does not do so at any uniform algebraic rate depending 
only on Eq, t, as may be seen by considering perturbations with support arbitrarily far from 
the shock location x = 0. See [Zlj for further discussion. 

It is a strength of this approach that such data may be treated nonetheless, and in a 
simple fashion parallel to the treatment of the constant-coefficient case. However, phase- 
asymptotic stability does not seem to be accessible by this simple argument scheme. For 
proofs of phase-asymptotic stability under strengthened assumptions on the initial data, 
involving additional pointwise information on the solution, see [Rf IHZ[ IHRl IHRZ[ IRZ| . 

4 The system case 

We have described the nonlinear stability argument of \Z1\ IMaZ2t IMaZ4j in the simple 
scalar setting of Burgers equation. We now discuss briefly how this carries over to the 
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case of general hyperbolic-parabolic systems, including Navier-Stokes equations of com- 
pressible gas dynamics and MHD. Namely, Remark 13.41 plus essentially the same argument 
described here gives nonlinear orbital stability of viscous shocks provided that they satisfy an 
Evans function (generalized spectral stability) assumption yielding the necessary pointwise 
bounds. The Evans condition is necessary for linearized stability as shown in |ZmiMaZ3] . It 
holds always for small-amplitude shocks, but may fail in general for large-amplitude shocks. 
In the large-amplitude case, it is readily checked numerically; in certain special limits, it 
may be checked analytically using asymptotic ODE and or singular perturbation theory. 
When the Evans condition fails, there are interesting implications for dynamics/bifurcation; 
see [ZSllZTllTZTllTMlTMlTZil lSS llB^ . 

See jAGJl IGZ[ IZIH IZSl IMaZ3j for discussion of the Evans function and its origins. For 
verification of the Evans condition for small-amplitude shocks, see |ZHl IHuZj \PZ\ |FS1| . For 
examples of unstable shocks, see [GZl IZSj . For numerical and analytical verification for 
large-amplitude shocks, see plZj IBHRZl iHLZl [HLyZj [CHNZ]; see [El EEZl iBDCTl [HiIZ2] 



for more general discussion of numerical Evans function techniques. See 
IGMWZl IGMWZ21 |FS2] for extensions to multiple dimensions. 

The derivation of pointwise Green function bounds for general systems is complicated, 
involving detailed estimates on the resolvent kernel using Evans function and asymptotic 
ODE techniques, converted to bounds on the Green kernel via stationary phase estimates 
in the inverse Laplace transform formula. See [ZHl ESI [Z2l EH [BeSZj K^MWZl l(^MWZ2j 
for discussions of these and related techniques. These are details of the linear theory. Here, 
we have chosen to isolate the nonlinear iteration argument by restricting to a case (Burgers 
equation) for which the linear theory is explicitly known a priori, in order to give the reader 
a flavor of the arguments. 

We emphasize: once the linearized theory is established, the nonlinear shock-tracking 
argument of \Z1\ \MaZ^ MaZ4 , is essentially the same for system or for scalar case. 
See Remark 13.41 

Acknowledgement. Thanks to Mark Williams and Benjamin Texier for their interest 
in the work, and for several helpful comments improving the exposition. 

APPENDICES 

A The small-amplitude limit 

It is instructive to consider the small-amplitude limit \u+ — U-\ — > 0. Consider now the 
family of stationary viscous shock solutions 

(A.l) u^(x) := — etanh(ex/2), lim = u'^ {x)u'^^ = 

x— >±oo 

of (jl.ip . and examine behavior as e ^ 0. 

Denote the associated homogeneous linearized equation by 

(A.2) vt - L^'v = vt + {a%x)v)x - Vxx = 0, v\t=o = f 
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where a^(x) := u'^{x). The invariance {x,t,u) — > {x / e,t/ ,u/ e) of Burgers equation con- 
verts this to the e- independent case (jl.Sp considered in Section [3l from which we may deduce 
the e-dependent Green function formula 

/+00 
G'{x,t-y)f{y)dy, 
-oo 

where 
(A.4) 

G^{x, t; y) := e^'%{x) = (n^)'(x) (1) (errfn( ^ ^ - errfn( ^ ^ 



+ ( I — -] + - 



> e 2 + e 2 / v 47rt ^e2+e 2/ v 47rt 
we are using the scahng relations 
u'^ (x) = eu{ex) and G^{x,t;y) = eG{ex,e'^t;ey). 
Decompose again 

(A.5) G%x, t; y) := E%x, t; y) + S%x, t; y) + W{x, t; y), 

where 
(A.6) 

E%x,t;y) := {u')' ix)e%y,t), e%y,t) := ( J-) (errfn(^^) - errfn( ^ + 



2eJ\ v/4t V4t 



(A.7) 5"(x,i;y):=^^ —) 7^^ + 



' e 2 + e 2 / V 47rt ^e2+e 2/ v 47rt 
and 



R%x,t;y) := - j (^errfn( J_ ) - errfn(^^) 



(A.8) 

+ — ) ( errfn( — —7= — ) — errfn( 



^2eyV ' Vit V4t 
Defining the perturbation 

(A.9) u{x,t) := u{x - a{t),t) - u'^{x), 

setting G^ := S"^ + i?^, and following the steps of Section [3l we obtain again the integral 
representation 



u{x,t)=l G{x,t;yyuo{y)dy 

(A.IO) 



00 

t roo 



y 

lO J -00 



Gl{x,t- s;y){N{u) + du){y, s)dy ds, 



15 



a{t) 



(A.ll) 



+ 



— oo 

t /•+00 



^-oo 



^Itiy^t - s){N{u) +au){y,s)dy ds. 



a{t) 



(A.12) 



+ 



— oo 

t /■ + 00 



e%y,t)uo{y)dy 
eliy, t - s){N{u) + a u){y, s)dyds. 



Dependence on e. Evidently, we could carry through the entire stability analysis of 
Section [3l as the e-dependent Green function = + + has the same form as 
G. However, the bounds obtained in this way- in particular, the estimate (j3.1ip on the 
remainder R- would involve constants C = C(e) > blowing up as e — > 0. This means 
that the allowable size £"0 < ^Q^^^y^ of perturbations, determined in the proof of Corollary 
12.41 goes to zero as e — > 0. That is, the basin of attraction of the shock established by our 
basic stability argument shrinks to zero as e ^ 0. Indeed, the bounds derived for general 
systems in [ZHl IMaZ3] (described briefly in Section [4]) share this same property, and so the 
basin of attraction for the stability results proved in [MaZ21 IMaZ4l \Z1\ \Z2\ IHZj and related 
works go to zero as the shock amplitude goes to zero. 

However, this is not an inherent limitation of the method, or the shock. Following, we 
show that by different, more careful, estimates of E'^ and i?^, we may in fact recover a 
uniform stability result, valid for perturbations of sufficiently small size independent of e. 



Proposition A.l. For some C > independent of e, < e < 1, and all t > 0, 

f +00 

G'{x,t;y)f{y)dy 



(A.13) 
(A.14) 



-00 
+00 



Gl{x,t;y)f{y)dy 



LP{x) 



LP{x) 



<ct-^(i-i/p)|/lii, 



<Ct-i(i-Vrt-||/|^,. 



+00 



&{x,t-y)f{y)dy 



LP{x) 



<C\f\ 



LP, 



+00 



Glix,t;y)fiy)dy 



LP{x) 



<Ct-2\f\LP; 



/+00 r+00 
eKy,t)fiy)dy < Cr^l/]^,, / ely{y,t)f{y)dy < Ct" VIl^ 

-00 J —00 



(A.16) 



+00 



el{x,t-y)f{y)dy < 



+00 



<Gt-2\f\L^- 



and 
(A.17) 



/+00 f+00 
e'{y,t)f{y)dy < Ce'VlLi, / el{y,t)f{y)dy < Ce-^-i/VlLi 

-00 J —00 
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Proof. As evidently obeys the same decay estimates as S, to establish the stated bounds 
on G^, it is sufficient to establish them for R^. This is a straightforward consequence of 
Lemmas IC.4I and IC.5I established in Appendix IC.3i Likewise, for the stated bounds on 
it is sufficient to establish corresponding U'{y) bounds on e^, from which the results then 
follow by Holder's inequality. The needed bounds are established in Lemma lC.31 Appendix 
[021 □ 

Corollary A. 2 (Stability of small-amplitude shock solutions). For < e < 1, viscous 
shock solutions u^{x) of (jl.ip are nonlinearly stable in n L°° and nonlinearly orbitally 
asympotically stable in 17 , p > 1, with respect to initial perturbations uq with L^r\L°° norm 
less than or equal to rjQ > sufficiently small, where rjo is independent o/O < e < 1. More 
precisely, for some C > independent o/O < e < 1, there is a ^ W^'°°{t) such that 

\u - u%- - a)\Lp{t) < C(l + 

(A.18) |d(t)| <C(l + t)-5i?o, 

\a{t)\ < Ce~^Eo, 
\u - M^liinL-'W < CEq, 
for all t > 0, 1 < p < oo, for solutions u of (jl.ip with Eq := \u — ^i'^|LinL°° \t=o < Vo- 

Proof. The proof of the first two bounds follows exactly as in the proof of Corollary 13.71 
in the fixed-amplitude case, since the integral equations for (n, q) form a closed system 
involving only G^, ef and e^j, and the bounds on G^, ef and e^j are the same as the bounds 
on on G, et and Cyt in the fixed-amplitude case. With these bounds established, we obtain 
the third bound from (IA.12|) . using the fact that the bounds on and are no worse than 
times the bounds on e and Cy in the fixed-amplitude case. 

Finally, we note that — u'^{x) = u{x — a{t),t) + {u'^{x) — uf{x — a(t)), so that 

\u{-,t) — is controlled by the sum of \u\ and \u^{x) — 0{x — By monotonicity of 

scalar shock profiles as orbits of the first-order scalar profile ODE (II. 4p . u^{x) — u'^ {x — a{t)) 
has one sign, hence 

/+00 
(■li^(x) — u^{x — a{t)))dx = \a{t)\\u\ — 
-oo 

and, by (fXTSjl (iii) . 

\u''{x) - - a(t))|ii = 2e\a{t)\ < 2GEq. 
Likewise, by the Mean Value Theorem, 

\u^{x) - Tf{x - a{t))\ < \a{t)\{u^)'\Loo < {GEo/e){e^) = CEqe, 

by the asymptotics uF)' ~ e^e~^^l^'L Thus, \u^{x) — u^{x — a{t))\ii^ioa < CEq, and so 
\u{x,t) — u'^{x)\x,ir\L°° — CEq for all t > 0, for Eq sufficiently small. This verifies the fourth 
inequality, yielding nonlinear stability and completing the result. □ 



17 



Remark A. 3. In the small-amplitude limit e — > 0, the shock shift a — > +00 as e times 
perturbation mass. Nonetheless, the stability estimates are uniform, independent of e. 

Remark A. 4. As discussed in Section 13.41 we have obtained stability for a class H 
of perturbations that lead to shock shifts a not only of order but also decaying 
subalgebraically to their limits a(+oo), if they exist. 

Remark A. 5. Here we have treated only the simple and explicit case of Burgers equation. 
It would be very interesting to try to treat the small-amplitude system case by a similarly 
simple argument based on this approach, using the singular perturbation techniques devel- 
oped in |MaZ3[ IPZ] to obtain the necessary sharpened e-dependent bounds analogous to 
those of Proposition lA.ll in the Burgers case to try to obtain results uniform in e. 

B Alternative shock-tracking schemes 

As discussed in Remarks 13.21 and 13. 4^ the quantity a{t) introduced for technical reasons in 
(|3.19p . has an interpretation as an "instantaneous shock shift", measuring the approximate 
location of a perturbed viscous shock profile at time t. This suggests the question what 
is the "exact" location of an asymptotic shock profile, and how well a{t) approximates 
this location. The study of this question leads to an interesting class of alternative shock- 
tracking schemes that are time-asymptotically equivalent to p.20p - (|3.2ip . based on localized 
projections, converging in the "infinite-localization" limit to a pointwise phase condition 
introduced in [GMWZj in the context of the small- viscosity limit. 

Unlike a perturbed inviscid shock wave, which is sharply located by the presence a 
discontinuity, a perturbed viscous shock wave is smooth, so requires some extrinsic criterion 
to define its location. Two intuitive definitions immediately come to mind. The first, 
defining the location of an unperturbed stationary scalar shock u = u{x) without loss of 
generality to be the origin, x = 0, is simply to define the location a{t) of a perturbed shock 
u as the point a{t) at which u takes on the value u{0), or 

(B.l) u{a{t),t) =u{0). 

By the Implicit Function Theorem and the fact that u'{0) 7^ (recall that u is monotone, 
as the solution of a scalar first-order traveling- wave ODE), this uniquely defines a for 
\u' — u'\L°°{t) sufficiently small. 

In the system case u G M", we cannot satisfy (jB.ip for all n coordinates using the single 
parameter a, so we must choose some preferred coordinate direction, substituting for (|B.ip 
the system analog 

(B.2) i-u{a{t),t) =£-u{0) 

for some vector i E such that i ■ u'{0) ^ 0, a condition that, by the Implicit Function 
Theorem, guarantees that a{t) is well-defined for \u' — n'|ioo(t) sufficiently small. 
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Defining the perturbation variable 



(B.3) u{x , t) = u{x + a{t) , t) — u{x) , 

following the notation of Section [3l we find that ()B.2p translates to the phase condition 
(B.4) £-M(0,t)=0, 

determining a{t) implicitly through ()B.3p . Condition ()B.4p is particularly natural from 
the point of view of the resolvent equation arising in solution by Laplace transform of the 
associated linearized equations. For, the resolvent equation consists of an underdetermined 
ordinary differential boundary- value problem for which the standard treatment is to remove 
indeterminacy by one or more phase conditions like ()B.4p . Indeed, this condition was 
introduced in [GMWZj starting from just such considerations, for the study of shock stability 
in the vanishing viscosity limit H 

The second intuitive definition is, following Goodman [G], to define the shock shift a so 
as to minimize the least squares distance of u{x, t) from the shifted shock u{x — a{t), that is, 
to minimize \u{-,t)\i^2. This leads to the "localized projection condition" (Euler-Lagrange 
equation) 

(B.5) {e,u)L2=o, {e,u')L2 = h 

where £{x) := (see Appendix IC.4I for this calculation). Here, the word "localized" 

I" 1^2 

refers to the fact that i{x) decays as a; ^ ±00. More generally, we denote as a localized 
projective condition any condition of form (IB. 51) with £ £ L^. This can be viewed as a 
nonlocal version of the pointwise phase condition ()B.4p . converging to ([Rl]) in the "infinite- 
localization limit" £{x) — > £o6{x), £q G M" constant, of a Dirac measure. 

Each of these schemes (either of form (1B.4P or (IB.SP ) may be written as an evolution 
equation in (u,a). Defining the perturbation variable u of ()B.3p . we find as in Section [3] 
that u obeys the partial differential equation 

(B.6) Ut — Lu = N{u)x + aiux + u^) 

depending on a, defined implicitly by (|B.5p . Differentiating (jB.SP with respect to t, we 
obtain 

= {£, ut)L2 = {£, Lu + N{u)x + a{ux + u^)) ^2 , 
which, using {£,Ux)l^ = 1, reduces to a{l + {£,Ux)L'2 = —{£, I-'U + N{u)x) l'^, or, rearranging, 

(n 7^ A, {£,Lu + N{u)x)l- 

1 + {£,Ux)l2 

well-defined for u G with \u\h2 sufficiently small. See [G llTZH IZ7j for related discussion. 



More precisely, a multi-dimensional version reducing to (|B.4|I in the one-dimensional case 
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Together, ()B.6P ~ ()B.7p determine a closed system of evolution equations for (n, d), similar 
in spirit to the system (j3.20p ~ (j3.2ip of Section [3l but local in time, whereas the system 
(|3.2Up - (|3.2ip involves "memory terms" depending on values of n, a at earlier times s < t. 
For each choice of test function i, there results a different evolution system, and different 
solutions (n, q) and a, representing different decompositions of the common solution u of 
(jl.ip under investigation, a perturbed viscous shock wave. 

We know already from the analysis of Corollary 13 . 71 that the solution u exists for all time, 
and converges to the set of translates of the background shock u. However, it is not a priori 
clear that the system (IB.6p -( lB^ has a global solution for any particular choice of i, nor 
that the solution u should decay as i — > 0. That is, it is not clear which of these alternative 
shock tracking schemes gives an accurate estimate of shock location in the sense that the 
known convergence of u to the set of translates is revealed by decay at the appropriate rate 
of the perturbation variable u. 

The following proposition asserts that all of these schemes are accurate in this sense, so 
that in principle any one of them could be used as the basis of an argument for nonlinear 
stability. Indeed, all lead to the same rates of decay. 

Proposition B.l. Let u'^'^^ , a^'^^ denote the solution of (j3.20p - (j3.2ip of Section \^ with 
initial data uq — u, Eq := |uo — 'ul^.inH^ sufficiently small, andu, a denote the solution with 
same initial data of ()B.6p - ()B.7p . with I ^ . Then, u, u^°^ exist for all t>Q, with 

\u\Lv{t), W'\Lv{t) < CEo{l + tr'2i^-'/P\ 

\u\H2{t), \u'''\H2{t)<CEo{l + t)-'/\ 
(B_8) Hi^nH^t) - \u'^'\Lir^H<t) < CEo{l + t)-^/^ 

\u - u\Lir,H^{t) < CEo, 
\am, la^^'mKCEo, 

\a-a'^^\{t) < CEo{l + t)-'^/^. 

Proof. A routine extension of the proof of Corollary 13.71 using the additional assumption of 
smallness of the initial data yields (j3.27p augmented with \u'^^^}{2{t) < CEq{1 + t)'^, 

We omit the details. (But see the results of [MaZ2t IMaZ4] in the much more complicated 

system case.) The corresponding bounds ()B.8P (i)-(ii). hence global existence of u, thus 

follow provided that we can estabhsh (lR8]) f 111 . 
Expanding 

u{x, t) = u{x + a{t),t) — u{x) 
(B.9) = u{x + a{t),t) - u{x + (a - a'^^'O) + ^ix + {a - a''^'')) - u{x) 

= u'^^ix + {a- a'^^),t) + (n(x + (a - a'""^)) - u{x)), 



we find using the Triangle inequality, followed by the Mean Value Theorem together with 
exponential decay of u' , that 

Hi^nH^it) - W'^^lL^nH^t) < \u{x + {a- a'^'O) - u{x)\^,^^^ < C\a - a'^^\{t), 
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so that (jB.Sp fiii) follows from (jB.Sp fvi). Likewise, (iv) follows from (i)-(iii) and (v), which 
in turn follows from (vi) and the bounds on |a''*'^|(t) established in Section [3l 

Thus, it remains only to prove (jB.Sp fvi). Applying definition {£,u)i2 = to expansion 
()B.9p . we obtain 

{^y^'ix + (a - a-^),t))^2 = + (a - a^^)) - u{x)) 

= -{l,{a-a'''^))v! + 0{\a-Q'^^)\^))L2. 

Applying now {i,v!)i2 = 1, and rearranging, we obtain 

W - c^'^^'m < \e\L^{\u'''\L-{t) + C\a- a-T) 

yielding (jB.Sp fvi) provided \a — a^'^^l is sufficiently small. The result then follows by conti- 
nuity of a, a^^^ and smallness of a at t = for Eq small, recalling that a^'^^(O) = 0. □ 

Remark B.2. As the only bound used on £ was its norm, the proof of Proposition lB.il 
is easily adapted to the case that £ is a bounded measure, in particular the case of a phase 
condition (IB.4p . This includes also more general cases such as the sum of point measures, 
leading to a sort of "difference stencil" condition determining shock location. 

Remark B.3. Recalling that a^^^{t) in general decays at most at subalgebraic rate (see 
Remark 13. 9p . we see from ()B.8P (iv) that a and a^'^^ are time- asymptotically equivalent in 
the sense that \a — a^^^\ decays at a rate faster than the (general) rate of decay of ja'^^^l. 

Remark B.4. For initial data in addition decaying as |uo(2;)| < CEo{l + it is 

shown for general systems in \iiR\ IRZ| that a^'^^ decays at the faster rate 

(B.12) la'^'^it) - a''^^(+oo)| < CEq{1 + t)-^/^. 

However, the same analysis yields sharpened bounds on u^^^ as well, giving also 

\u'^^{x,t)\ <CEo{l + ty^ for \x\ <9t, 

9 > sufficiently small. Substituting in (jB.lOp . we obtain in place of ()B.lip the estimate 



a 



a^'^^Kt) < CEo{l + t)-^\i\Li + CEo{l + f \e{x)\dx + C\a - a'"^^^ 

J\x\>et 



yielding \a - a'^^it) < CEo{l + t)-^ provided \£{x)\ < C(l + |x|)-3/2. 

Thus, under this strengthened decay requirement on i, we obtain time-asymptotic 
equivalence of a and a^'^^ also in this case. Bound (|B.12p is sharp, as can be seen by direct 
computation on the linear term in (j3.19p for data decaying as (1 -|- (Note that 

the linear 0{Eo) term dominates the nonlinear 0{Eq) term up to any finite time, for Eq 
sufficiently small.) 
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Conclusions. By comparison with the scheme of Section [3l we find that each of the 
alternative shock-tracking schemes described in this Appendix, based on locahzed phase 
conditions, yields a globally defined solution exhibiting the same rates of decay as the 
perturbation w^*^^ defined in Section [3l That is, essentially any tracking scheme based on in- 
formation that is "local to the shock" in the sense that it is accessible by inner product with 
an function (resp. bounded measure) i yields a convergent system of perturbation equa- 
tions. Note, further, that the only information used to draw these conclusions consists of 
estimates on {u'''^,a"'^) already established in [ZTl IMaZ2l [MaZ3l IMaZ4l [HZ. RZ] for Evans- 
stable Lax or undercompressive type shocks of general hyperbolic-parabolic systems. Thus, 
the conclusions of Proposition iRJl and Remarks \B.SA\K4\ remain valid for Evans-stable Lax 
or undercompressive shocks of general systems of hyperbolic-parabolic conservation laws^ 

An interesting question is whether we could carry out a nonlinear stability analysis for 
these schemes from first principles rather than by comparison to our existing results. This is 
particularly intriguing for the case of the pointwise phase condition (|B.4|1 , for which resolvent 
(and thus pointwise Green function) bounds are available through the framework developed 
in |GMWZ| . Besides the intrinsic interest of this question, there are real advantages to the 
scheme based on ()B.4p for extension to more complicated situations: for example, the fact 
that it is local in time (the scheme in Section [3] by contrast involves "memory terms ") , and 
that the phase condition ()B.4p makes no reference to the explicit structure of the system. 

C Miscellaneous estimates 
C.l Bounds on e 

Lemma C.l. For some C > and all t > 0, 



(C.l) \e{;t)\L^,<C, 

(C.2) \ey{;t)\LP,\et{;t)\LP <Ct-'2i^~'/P\ 

(C.3) |ei,(-,t)|LP <Ct-5(i-Vrt-i/2^ 

(C.4) \ey{y,t)l\et{y,t)\<Ct-'^'[e-^^ +e-^^), 
(C.5) \ety{y,t)\<Ct-'[e-^^ +e-^^). 



^ With the inclusion of additional phase conditions to account for additional degrees of freedom in the 
time-asymptotic state (see |HZ. .RZ^ l. these methods and estimates extend also to the overcompressive case. 
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Proof. Bound (jC.ip follows immediately from definition (j3.7p . Given ()C.4p - ()C.5p . bounds 
()C.2p - ()C.3P follow as in the heat kernel estimates (j2.7p - (j2.8p . Thus, it remains only to estab- 

lish (|U:4l) - (|U3D . Differentiating we have ey{y,t) = i:;^^] ^ 

yielding (|C.4p (i). Differentiating (j3.7p with respect to t, we obtain 



e 4t 



et{y,t) 



(C.6) 



-1 



u+ — 



e it e 4t 
+ 



i-y - 1) e 



{-y-tf 



(-y + t) e 4t 



yielding (|C.4p fii) immediately for t > 1. By the Mean Value Theorem, for t < 1, 



(-y - 1) e 



(-y-tr 



{-y+tr 



{-y + t)e it 



Vt 



'A-Kt 



-1 



z e it \ 



z=-y+6t 



de 



(C.7) 



< 2Ct 



{-y-tr 



{-y+tr 



< Cie ct + e ct 



which, together with (jC.6p . yields again (|C.4p (ii). Estimate (jC.Sp goes similarly. Note that 
we have taken crucial account of cancellation in the small time estimates of et, ety. □ 



Remark C.2. For t < 1, a calculation analogous to (ICTTj) yields \ey{y,t)\ < Ce 
and thus |e(-, ^ as s — > 0. 

C.2 Bounds on 

Lemma C.3. For some C > 0, all < e < 1, and all t > 0, 



(y + a_t)^ 



(C.8) 



\e'i-,t)\Loo,<C/e, 



(C.9) 

(C.IO) 

(C.ll) 

(C.12) 



\e%{;t)\L.<{C/e)t-^2(^-'M, 



ty 



(•,t)|L. <Ct-5(i-Vrt-i/2^ 



\e%{y,t)\<{C/e)t-'/' 



e ct ^ e Ct 
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(C.13) 



(C.14) 



\e'tiy,t)\<Ct-'/^e-^^ +e 



(-y+tr 

Ct 



\e'tyiy, t)\ < Cr^ e — + e — 



Proof. Bounds (|C.8h . (IC.9p . and (|C.12p follow exactly as in the e-independent case. Bound 
(jC.lOP follows immediately provided that we can establish ([aT3]l . as we now do. Differen- 
tiating (jA.6P with respect to t, we obtain 



e?(y,t) 



(C.15) 



-1\ / e 4t e 4t 
ft ' \ fi-y - t) e 4t — (_y + et)e 4i— 



2e ^ V ^ ^/i^ ^ ^/i^ 

yielding (|C.13P immediately for t > e^^. By the Mean Value Theorem, for t < 



(— y — rt) e 4t (— y + et) e it 



Vt 



Vt 



(C.16) 



< 2Cet 



1 



2=— y+0et 



5. 



z e 4t 



\z=~y 



<Ce[e Ct + e ct 



which, together with (|C.15[) . yields again (IC.lSp . Bounds (|C.11|) and (|C.14|) follow similarly. 



□ 



C.3 Bounds on W 

Lemma C.4. For fCf := K{x, y)f{y) dy and any 1 < p < oo, 



(C.17) 



\ICf\Lp < sup \K{-,y)\Lp\ fill, 
y 



(C.18) 

Proof. By the Triangle inequality 



\}C\lp^lp < max{su-p\K{x,-)\Li,sup\K{-,y)\]^i} 

X y 



K{;y)f{y)dy < \K{;y)\Lp\f{y)\dy < sup 
LPi^) Jr y 



establishing (jC.17p . This yields also (jC.lSp in case p = 1. Likewise, 

K{x,y)f{y)dy < / y)|d?/|/|Loo < sup OIli 



establishing the (|C.18p for p = oo. For general p, (|C.18p then follows by the Riesz-Thorin 
Interpolation Theorem. □ 
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Lemma C.5. For some C > 0, allQ < e <1, and all t > 0, 

(C.19) sup|i?^(-,i;y)|iP(,), sup t; OIlpc^) < Ct-^a-Vp)^ 



(C.20) 



Proof. From (u^)' ~ e^e ^'^'^l, we obtain 



(C.21) 



R'{x,t-y) = {l/2e)x{if)'{x) 



0x-y-£t)^ 
it 



e 4t 



(Bx-y + ety 

e 4t 



/o ^ v4vrt VAnt 
from which we obtain immediately |i?^|ioo < Ct^^^"^, and, bounding Ce~^^^^^ by C, 

sup|i?^|iP(^) <Ct-^(i-i/p) 

X 

for any p. 

Bounding the integral on the righthand side by Cit~^^'^ and the L^{x) norm of Ce"^^'^' 
by C2/e, we find sup^ {R'^Il^Ix) < C2t~'^^'^/e < C for t > For t < e~^, on the other 

hand, we may estimate the integral (the middle displayed term in the first equality) instead, 
using the Mean Value Theorem, as 



Bx-y-ety 



(Sx-y + sty , (Ox-y-zy 

e 4t \ r /e 4t 

^^]de< / i2et)d/ 

VAnt ) Jo 



\z=Zt.&[~et,et] 



de 



< {2et) I Cr^ d9 < Ce, 
^0 



to again obtain sup^ < C2£/s < C. The bounds on supj^ \R'^\lp{x) then follow by 

Holder interpolation between the and L°° bounds, verifying (jC.lOp Similar computations 
yield dU^QD- □ 

C.4 Euler— Lagrange equations for least squares 

Setting E{a) := ^\u\'j^2 = f |^^(" + o.,t) — n(-)|^2 and differentiating, we have 



dE 
da 



{u{- + a,t) - u{-),u\- + a,t))L2 = {u,u' + u')l2 = {u, u')l2, 
where, in the final equality, we have used {u,u')i^2 = j^^{v? /2)' {x)dx = for u G . 
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